Fluctuations of a weakly interacting Bose-Einstein condensate. 



00 
O 

o 

(N 
o 

Q 

m 



in 



i 

C 

O 

o 



> 
(N 

O 

q 

00 

o 
i> 
o 



X 



Zbigniew Idziaszek, 1, 2 Lukasz Zawitkowski, 2 Mariusz Gajda, 3,4 and Kazimierz Rzazewski 2,4 

1 Institute of Theoretical Physics, University of Warsaw, 00-681 Warsaw, Poland 
Center for Theoretical Physics, Polish Academy of Sciences, Aleja Lotnikow 32/46, 02-668 Warsaw, Poland 
:l Institute of Physics, Polish Academy of Sciences, Aleja Lotnikow 32/46, 02-668 Warsaw, Poland 
^Faculty of Mathematics and Sciences, Cardinal Stefan Wyszynski University, Warsaw, Poland 

Fluctuations of the number of condensed atoms in a finite-size, weakly interacting Bose gas 
confined in a box potential are investigated for temperatures up to the critical region. The canonical 
partition functions are evaluated using a recursive scheme for smaller systems, and a saddle-point 
approximation for larger samples, that allows to treat realistic size systems containing up to N ~ 10 5 
particles. We point out the importance of particle-number constrain and interactions between out 
of condensate atoms for the statistics near the critical region. For sufficiently large systems the 
crossover from the anomalous to normal scaling of the fluctuations is observed. The excitations are 
described in a self-consistent way within the Bogoliubov-Popov approximation, and the interactions 
between thermal atoms are described by means of the Hartree-Fock method. 
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A breakdown of the standard, grand-canonical ensem- 
ble to describe fluctuations of an ideal Bose gas and a 
necessity for canonical or microcanonical description has 
been noticed already long time ago [l[ , but only in recent 
decade the problem of fluctuations has received renewed 
attention due to the experimental achievement of Bose- 
Einstein condensate (BEC) in ultracold trapped gases. 
For ideal gases, the canonical and microcanonical fluctua- 
tions have been thoroughly investigated 
SUB], and severalpowerful techniques, like the Maxwell 
Demon ensemble (j,[6|,LZ|: have been developed. For inter- 
acting particles the fluctuations have been studied mainly 
within the Bogoliubov approximation ll| of weakly in- 
teracting gases 0EIIll0E [ljtlrat proved to 
be extremely successful to describe many other proper- 
ties of BEC. The exact treatments, so far applied only 
for one-dimensional systems [la ], confirmed an excellent 
agreement with predictions of the Bogoliubov method. 
We note that some controversy exists about the applica- 
bility of the mean- field theory to this problem [l|| , on the 
other hand, other approaches, like the perturbation the- 
ory, lead to qualitatively different results for fluctuations 
of relatively small condensates 2^, 21 1 . The ultimate ver- 
ification will be done in experiments. However, to date 
only the statistics of the total number of atoms has been 
measured [22], and a technique involving scattering of 
short laser pulses has been proposed [23| but not realized 

So far the studies of fluctuations in weakly interacting 
gases have been limited to the regime of low tempera- 
tures, and only recently the critical region (close to the 
critical temperature T c ) in a finite-size system has been 
explored 24j. In this case the Bogoliubov-Popov (B-P) 
approximation [25| has been applied to account for the 
condensate depletion at finite temperatures and to obtain 
a description that smoothly interpolates between the de- 
generate regime below T c and an ideal gas statistics above 



In this Letter we reinvestigate the problem of fluctu- 
ations for weakly interacting gas, putting a special em- 
phasis on the interactions of out of condensate atoms, 
that apart from the critical region, turn out to be im- 
portant even at moderate temperatures. Following the 
Bogoliubov-Popov approximation for a uniform Bose gas 
of N atoms confined in a three-dimensional box of size 
L with periodic boundary conditions we start with the 
Hamiltonian: 



H = H F 



E ex (N,N Q ) 



E 



ek&k&k + E e 



(N,N ). (1) 



Operators 6k = U^a^ + Vka^ k are the Bogoliubov quasi- 
particle annihilation operators, obeying Bose commuta- 
tion relations [6k, 6 k ,] — $k,k', &k represent annihilation 
operators for a mode with quantized momentum tik. The 
celebrated Bogoliubov-Popov energy spectrum: 



£k = \J (e k + gn Q ) 2 - (gn y 



(2) 



depends on the condensate density no = No/V. Here, 
e k = 4n 2 h 2 k 2 /mL 2 is the kinetic energy of a mode k, 
to is the mass of atoms, g — Ai:h 2 a/m is the interaction 
strength, and a is the s-wave scattering length charac- 
terizing the contact potential V(r — r') = gS^^r — r'). 
Bogoliubov coefficients satisfy equations: + Vj 2 = 
(gno + e k )Ak = W k and U 2 - V 2 = 1 [25]. Finally, 
E ex (N, No) describes the interaction energy between out 
of condensate atoms, which in the B-P model can be cal- 
culated on the level of Hartree-Fock (HF) approximation: 
E ex (N,No) = ^N 2 ex [H, with N ex = N - N denoting 
the number of atoms in excited (k 7^ 0) modes. This 
corresponds to taking only secular part of interactions 
between thermal atoms. The considered Hamiltonian ne- 
glects a finite life-time of quasiparticle excitations arising 
from interaction between quasiparticles [2^ ]. 
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The canonical-ensemble partition function for a system 
with N atoms and temperature ksT = is 



TV 



z( N ,p)= £ £...£•• 



(3) 



7V ex =0 ni=0 



where nk are populations of quasiparticle excitations, 
E = Ek#o £ k(iVo) n k + E ex (N ex ) is the energy and 

N ex = Sk^o n kWk + ^k is the number of thermal atoms 
of a given configuration of excitations, that differs from 
the total number of excitations due to the Bogoliubov 
transformation [32| ■ Although the condensate mode does 
not appear in the sum in the Hamiltonian, its population 
affects the energy spectrum and the interaction energy of 
atoms in excited modes. In order to enforce the constrain 
on the total number of particles rigorously, we keep the 
energy spectrum dependent on the actual number of con- 
densed atoms, as follows from Eq. ©■ We calculate the 
conditional statistical partition function: 



Z No (N ex )= ••• E ■■•e- /3 ^ 



N ex ,N ex > 



(4) 



711=0 



"k=0 



which corresponds to a case with No condensed atoms 
and N ex thermal atoms. In terms of these functions the 
probability of finding N condensed atoms is P(N ) = 

Z ^ N - No) and Z = £^ o=0 Z No (N N ). 

The recurrence algorithm used in our calculations is 
an enhanced version of the earlier algorithm applied to 
the ideal Bose gas (IBG) [6], and it will be presented in 
details elsewhere. It makes use of the fact that Zm (N &x ) 
treats the number of condensed and thermal atoms as in- 
dependent variables and the number of condensed atoms 
becomes a parameter. As an intermediate step one ob- 
tains the following result for the mean number of quasi- 
particle excitations in mode q provided that there are Nq 
condensed atoms and N ex thermal atoms in the system 



I N 



i=i 



_ 0Ua Z No (N ex - lW k ) 
Zn (N ex ) 



(5) 



We calculate the mean condensate population (No) = 
J2n =o P( n o)N„ and its fluctuations (S 2 N ) = (N^) - 
(No) 2 . In the canonical ensemble (Nq) = N — (N ex ) 
and (S 2 N Q ) = (S 2 N ex ). The fluctuations, can be written 
as a sum of two contributions: (5 2 N ex ) — (S 2 N ex )t + 
(5 2 N ex )Q. The first term represents thermal fluctuations, 
that we calculate from the probability distribution P(Nq) 

(N 2 ex )r = £ W * W i «n k n q ) - (n k )(n q » , (6) 

N 

= ^2 N oP(No)~(N ). (7) 

7V =0 



The second term, (S 2 N cx )q describes the quantum part 
of the fluctuations, a non vanishing component at T = 
in an interacting gas [l2j]. They result from the Bogoli- 
ubov transformation applied to the quantum average of 
the N 2 X operator: 

m 2 ex ) Q =4£([/X 2 (n k n_ k +n k +i)) - ( 8 ) 

In the equations above, the average of an arbitrary op- 
erator can be expressed in terms of conditional aver- 
ages: (X) = En =o(X)n„~ N ° P ( N o), with the mean oc- 
cupation numbers (nq)^^ given by (0, and the correla- 
tion of modes with the opposite momenta (rikfi-k)]^" = 
Ey=i e-M+^«Z No (N ex -(1+ j)Wq) /Z No (N ex ). ° 

From the practical point of view the recursive method 
is applicable for systems of maximum a few hundred par- 
ticles. For larger TV, the calculations become numerically 
very demanding, and to treat larger samples we have 
developed a semi-analytical approach, that is based on 
saddle-point approximation to the contour-integral rep- 
resentation of Z(N ex ), known in the literature as Darwin- 
Fowler method 28] . Derivation proceeds basically in the 
same manner as for an ideal gas [J, |8( and yields 



Z No (N ex ,(3) 



2/Vo(zo,/3) 



/ d 2 



(9) 



where Sjv (z, 0) is the grand-canonical partition function 
for the excited subsystem, zq — e A ° denotes the posi- 
tion of the saddle point, determined by (N ex ) G c — N ex , 
with (N ex ) GC = -d\ \nE No (z , (3) denoting the grand- 
canonical expectation value for the number of excited 
atoms. In analogy to the ideal gas, Stv (z,/3) can be 
written in a closed form 



Zn (z,I3) 



n 



[1 



z Wk exp 



R3e k )] \ (10) 



A similar saddle-point method can be applied to deter- 
mine (rik) and (nk«--k) entering formula ([5]) for (5 2 Nq). 

This way we have obtained a scheme that allows us to 
calculate statistical properties of the weakly-interacting 
condensate at all temperatures. While we keep only 
the HF contribution to interactions between quasiparti- 
cles we otherwise preserve the number of atoms through- 
out the calculations, which requires inclusion of the en- 
ergy spectrum dependent on the actual number of con- 
densed atoms. This can be contrasted to the common 
approximation assuming the energy spectrum dependent 
on the mean number of condensed atoms: ek(Ao) = 
e k((A^o)) [lil HH Hil ' use d m a simple but useful models 
of thermal equilibrium with Bose-populated excitations 
[lil I20I] . This way one would obtain different formulas for 
Z = En =o Z(n o) (N - No) and P(N ) = W^ZfM, 
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FIG. 1: (Colors on-line) Logarithm of the canonical par- 
tition function as a function of No and N ex (a), along with 
cross-sections yielding the probability distributions of the con- 
densate (b) for rigorous (blue) and average spectrum (red). 
Parameters are: N = 1000, an 1/3 = 0.03 and T = 0.8T C . 

which, being numerically less demanding, require self- 
consistent determination of (No). However, such a seem- 
ingly natural simplification leads to observable distor- 
tions of the results. This is illustrated in Fig. [T] pre- 
senting the canonical partition function of the excited 
subsystem in the parameter space (No,N ex ). Inclusion 
of the excitation spectrum dependent on the actual num- 
ber of condensed atoms, correspond to performing a cut 
along a line Nq + N ex — N, whereas the approximation 
assuming average spectrum corresponds to a cut along 
No = (N ) = const, with (Nq) determined in a self con- 
sistent way. These two approaches yield probability dis- 
tributions of the number of condensed atoms (see right 
panel) that differ both in the position of the maximum 
and the width of the peak, that determine the values of 
(Nq) and (S 2 N ), respectively. 

The comparison of the mean condensate populations 
and fluctuations, for relatively small system of N = 200, 
and an 1 / 3, — 0.1, calculated for rigorous and average 
spectrum is presented in Fig [2] The results obtained 
in the model with rigorous spectrum and the HF con- 
tribution to a thermal atoms interaction differ substan- 
tially from the mean occupation and fluctuations ob- 
tained when the thermal atoms interaction is totally ne- 
glected. On the other hand we note that the two other 
approaches for an interacting gas lead to rather similar 
results, and some discrepancies can be only observed in 
behavior of fluctuations close to the critical temperature. 

The situation changes, however, for larger systems (see 
FiglHJl. For sufficiently large system of N = 10000 atoms, 
and an 1 / 3 = 0.05 inclusion of the rigorous spectrum to- 
gether with E ex significantly affect the condensate statis- 
tics. This is more evident in the case of fluctuations, that 
in our model remain much smaller than the fluctuations 
calculated in the model assuming average spectrum, even 
at temperatures much smaller than the critical one. 

Finally we have verified how the mean condensate pop- 
ulation and its fluctuations depend on the size of the sys- 




T/T c 

FIG. 2: (Colors on-line) Condensate population and fluctu- 
ations (inset) versus scaled temperature T/T c for a system 
of N = 200, and cm 1 ' 3 = 0.1, described by the Bogoliubov- 
Popov Hamiltonian Hb with (black dot-dash) and without 
(red dots) inclusion of interactions between of out of conden- 
sate atoms E ex . Blue solid line represents a corresponding 
ideal gas, while green dash line shows the results for the model 
assuming average spectrum Sk({iVo)) without E ex term. 
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FIG. 3: (Colors on-line) Condensate population (upper 
panel) and fluctuations (lower panel) versus scaled temper- 
ature T/T c for a system with N = 10000, an 1/3 = 0.05 ob- 
tained for: (black solid) Hb + E ex with iVo-dependent spec- 
trum, (red dash) Hb with average spectrum approximation, 
(blue thin solid) an ideal gas. 

tem, while keeping the interaction parameter an 1 / 3 fixed. 
Fig. [4] presents the results for the rigorous spectrum in- 
cluding the excited atom interactions, for art 1 / 3 = 0.05 
and the number of particles varying from N = 100 to 
N = 10 s . While the size of the system increases, fluctu- 
ations tend to be proportional to V^V, so they become 
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FIG. 4: (Colors on-line) Normalized condensate fluctuations 
{SNq)/N 1 ^ 2 and condensate fraction (inset) versus scaled tem- 
perature T/T c for fixed an 1 / 3 = 0.05 but varying N: (a, red 
dots) N — 100, (b, blue dash) N — 1000, (c, green dot-dash) 
N = 10 4 , (d, black solid) N = 10 5 . 



normal. On the other hand, for small systems the scal- 
ing remains anomalous. This result shows that anoma- 
lous scaling (5N ~ TV 2 / 3 ), predicted within the Bo- 
goliubov method neglecting the interactions of thermal 
atoms 12, [3 [3 12 1 1 holds only for relatively small 



number of atoms. 

One observes that for large systems the fluctuations 
exhibit a high and narrow peak close to T c . We are not 
sure of its physical significance since the B-P spectrum 
is questionable so close to the critical temperature. We 
point out that position of this peak, that define some 
characteristic temperature for our model, remain fixed 
in the thermodynamic limit (art 1 / 3 = const, N — > oo). 
We have verified that the position of the peak vary with 
interactions as AT c h/T c — l.bGan 1 ^ 3 , which is similar 
to the shift of the critical temperature in an interacting 
gas: AT C /T C = 1.29cm 1 / 3 [p. 

In this Letter we have presented the most complete 
to date discussion of the statistical properties of BEC 
confined to a box. In particular we have pointed out 
the importance of the strict enforcement of the parti- 
cle number conservation and of the interactions between 
thermal atoms. These two new elements, that have been 
neglected in the previous approaches, turn out to be rel- 
evant not only close to the critical region, but also at 
moderate temperatures, affecting the scaling properties 
of fluctuations in sufficiently large systems. 
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